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Abstract: In this paper, an amperometric glucose biosensor is modeled numerically. The 
model is based on non-stationary reaction-diffusion type equations. The model consists of 
four layers. An enzyme layer lies directly on a working electrode surface. The enzyme layer 
is attached to an electrode by a polyvinyl alcohol (PVA) coated terylene membrane. This 
membrane is modeled as a PVA layer and a terylene layer, which have different diffusivities. 
The fourth layer of the model is the diffusion layer, which is modeled using the Nernst 
approach. The system of partial differential equations is solved numerically using the finite 
difference technique. The operation of the biosensor was analyzed computationally with 
special emphasis on the biosensor response sensitivity to oxygen when the experiment was 
carried out in aerobic conditions. Particularly, numerical experiments show that the overall 
biosensor response sensitivity to oxygen is insignificant. The simulation results qualitatively 
explain and confirm the experimentally observed biosensor behavior. 

Keywords: modeling; reaction-diffusion; biosensor; oxygen; aerobic; anaerobic 



1. Introduction 



A chemical sensor is a device that transforms chemical information, ranging from the concentration of 
a specific sample component to total composition analysis, into an analytically useful signal. Chemical 
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sensors usually contain two basic components connected in series: a chemical (molecular) recognition 
system (receptor) and a physicochemical transducer. Biosensors are chemical sensors in which the 
recognition system utilizes a biochemical mechanism [1-3]. 

A biosensor recognition system is usually enzyme-based. Enzymes are known for their characteristic 
trait to catalyze only one specific chemical reaction, which provides biosensors with high specificity. 
Furthermore, enzymes are known to be very efficient catalysts, ensuring high biosensor sensitivity [4-8]. 

A physicochemical transducer may be developed by employing several different approaches. One of 
the most common approaches is amperometric. In this case, one of the products of chemical reactions 
occurring during the operation of a biosensor engages in an electrochemical reaction on the biosensor 
electrode surface. The measured current is usually proportional to the concentration of the analyte, which 
allows the determination of the analyte concentration using a preestablished calibration curve [2,5,9]. 

Carbon nanotubes are a popular electrode material among biosensor scientists. A lot of biosensors 
have been developed using carbon nanotube electrodes [10,11]. Some biosensors based on the 
direct electron transfer between an enzyme and a carbon nanotube electrode were designed [12,13]. 
However, it was demonstrated that it is possible to achieve direct electron transfer using a cheaper 
electrode material, modified carbon black [14,15]. In [15], a glucose biosensor employing the 
s-PQQ-glucose dehydrogenase (GDH) enzyme and the electrode based on carbon black was created. 
In this work, the glucose biosensor based on the same electrode and the same enzyme employing 
mediator N-methylphenazonium methyl sulfate was developed. To increase the stability and prolong 
the calibration curve of the biosensor, the enzyme was immobilized onto a polyvinyl alcohol (PVA) 
coated terylene film. Thus, the biosensor has the shape of a sandwich, where the permeable membrane 
entraps the enzyme layer [3,16] 

The understanding of intrinsic biosensor mechanisms is important in designing and optimizing novel 
biosensors. In order to fully understand the processes taking place during biosensor operation, a 
mathematical model of a biosensor should be built [17,18]. Various mathematical models of biosensors 
have been built and successfully used so far [19-21]. Different sandwich-type (multilayer) biosensors 
have been also mathematically modeled [22-25]. Comprehensive reviews on the mathematical modeling 
of amperometric biosensors have been presented [26,27]. 

PQQ-dependent enzymes do not react with molecular oxygen [15]; thus, the biosensors presented 
in [15] do not require anaerobic conditions during operation. However, the biosensor presented in this 
paper uses a mediator, which does react with molecular oxygen [28]. The goal of this paper is to assess 
the extent of oxygen's influence on biosensor operation if it is used in aerobic conditions. A mathematical 
model of a glucose biosensor presented in this paper has been developed recently [29] . The model did 
not consider the oxidation of a mediator by molecular oxygen present in the bulk solution. The new 
model was created in order to model the influence of oxygen on the biosensor response. 

The biosensor behavior was numerically analyzed at various values of input parameters of the model. 
The influence of the diffusion, as well as of the mediator's oxidation by oxygen on the biosensor response 
were thoroughly investigated. 
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2. Experimental 

Aiming to design a biosensor electrode powder, carbon black RAVEN-Mobtained from Columbian 
Chemicals Co. (Atlanta, GA, USA) was mixed with a pasting liquid consisting of 10% polyvinyl 
dichloride in acetone and further was extruded, forming a tablet [30]. The tablet was sealed in a 
Teflon tube. The electrode was washed with bidistilled water and dried before use. As a biological 
recognition element, soluble PQQ-dependent glucose dehydrogenase (sPQQ-GDH) from Acinetobacter 
calcoaceticus, E.C.I. 1.5. 2 was used. The sPQQ-GDH was isolated and purified by the method reported 
in [31]. The enzyme was immobilized on individual flexible supports of 0.1% polyvinyl alcohol 
coated terylene. 

The thickness of the terylene membrane was of 12 |am. A thin layer of the PVA was formed on the 
terylene membrane. It was estimated that the thickness of this layer was about 1 |am. 

All electrochemical measurements were performed using the electrochemical analyzer, 
PARSTAT 2273 (Princeton Applied Research, US), with a conventional three-electrode system 
containing the carbon paste electrode as a working electrode, a platinum wire as a counter electrode 
and an Ag/AgCl in saturated KCl as a reference electrode (all potential values presented in this paper 
are versus this reference electrode). The measurements were performed in potentio static conditions at 
E = 0.4 V. Acetate buffer (50mmol/L, pH = 6.0) was used as a default buffer. All measurements were 
carried out at an ambient room temperature (20 °C). 

The initial experiments were conducted in both anaerobic and aerobic conditions. However, the 
difference in the signal between anaerobic and aerobic conditions was not observed. Thus, the rest 
of the experiments were conducted in aerobic conditions. The data used in this paper are for experiments 
conducted in aerobic conditions. 

3. Mathematical Model 

3.1. Reaction Scheme 

We consider that the following chemical reactions take place during the operation of the 
biosensor [15,28,32,33]: 

h 

GDHox + glucose »-GDHrcd + gluconolactone (1) 

k2 

GDHred + PMSox ^GDHox + PMS^ed (2) 

k3 

PMSred + O2 ^PMSox + HO^ (3) 

PMS,ed ^ PMSox + 2e- (4) 

During the first chemical reaction, glucose dehydrogenase oxidizes glucose to gluconolactone. During 
the second chemical reaction, the reduced form of glucose dehydrogenase (GDHred) is oxidized by the 
mediator, N-methylphenazonium methyl sulfate (PMS), and regains its primary oxidized form (GDHqx)- 
The third reaction is the oxidation reaction of the mediator by the oxygen that is present in the solution. 
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During this reaction, the mediator is oxidized and regains its primary oxidized form. The fourth reaction 
is an electrochemical reaction that takes place on an electrode surface. During this reaction, the mediator 
is oxidized in the same way as in the third reaction. 

Reactions (3) and (4) are competitive, as they both are dependent on the same reactant, PMSrcd- 
A high rate of Reaction (3) may reduce the concentration of PMSrcd and, consequently, the rate of 
Reaction (4) and, thus, the electric current, which is the biosensor response. 

For the sake of simplicity, further in this paper, we use an abstract notation of chemical species. 
As the purpose of the biosensor is the measurement of the glucose concentration, glucose is called the 
substrate and denoted as S; gluconolactone is called the product and denoted as Pi; Eqx denotes GDHqx; 
Ered denotes GDHred; Mqx is PMSox; and M^ed is PMSred- P2 denotes the product of the third reaction: 
HO2 . Thus, the reaction schemes (l)-(4) transforms to: 

Eox + S ^Ered + Pi (5) 

Ered + Mox— ^Eox + Mred (6) 
Mred + O2 ^Mox + P2 (7) 

Mred ^Mox + 2e- (8) 



3.2. Biosensor Principal Structure 

The biosensor consists of three layers of different diffusivity of the species. The mathematical model 
should consider all these layers plus a diffusion layer, where concentrations of the substances differ from 
the ones in a bulk solution. In our mathematical model, we consider the Nemst model of a diffusion 
layer, which suggests that the diffusion front is stopped by the convection at a certain distance from 
the electrode. The profiles of concentrations inside a diffusion layer acquire linear shapes at a steady 
state. On the contrary, the semi-infinite model of the diffusion layer considers that the diffusion front 
may infinitely shift to the bulk of the solution. However, if the measurement time is not very short, it is 
indispensable to take into consideration the consequences of convection, as well [34-36]. 

Let us define di, d2, dz and ^4 as the thicknesses of the enzyme, PVA, terylene membrane and 
diffusion layers, respectively. We will also need values representing the distances between the electrode 
surface and the boundaries of the regions. Let ai, 02, 03 and 04 be the distances between the electrode 
surface and one of those boundaries, as shown in Figure 1. The thicknesses of all layers are specified in 
Table 1 . Thicknesses di and 6/4 were not measured during the experiment. Their values were estimated 
by the model. 

In the enzyme layer, biochemical and chemical reactions, as well as mass transport by diffusion take 
place. It is considered that the molecules of the enzyme are large enough and immobile; thus, diffusion 
does not influence the concentrations of both enzyme forms. Only the mass transport by diffusion and 
a chemical reaction (7) take place in the PVA layer, terylene membrane and diffusion layer, as enzyme 
molecules are absent in these three layers of the biosensor. 
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Figure 1. The principal structure of a biosensor. PVA, polyvinyl alcohol. 
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Table 1. The parameters of the model. 



Parameter 


Value 


Reference 




0 X lU 111 


eSlllllallOIl 


d2 


1 X IQ-^m 


experiment 


ds 


1.2 X 10"^ m 


experiment 


d^ 


1.5 X 10-^m 


estimation 


Ds,i, -Dmo,,i, -DM,,d,i 


1.5 X 10-1° m^/s 


estimation 


-Ds,2, -Dmox,2, -DM,cd,2 


4.2 X 10-1° m^/s 


estimation 


Ds,3, -Dmox,3; -DMrcd,3 


3.75 X 10-1° m^/s 


estimation 


DsA 


6.77 X 10-1° m^/s 


[37] 


-Dmox,4, -DM,.ed,4 


4.57 X 10-1° m^/s 


[38] 


-Do2,l) -Do2,2, -Do2,3 


1.970 X 10-9 m^/s 


estimation 


D02A 


1.970 X 10-9 m^/s 


[39] 


Co 


1 X 10-^mol/m^ 


estimation 


mo 


5 X 10-^ mol/m^ 


experiment 


So 


4.98 moVm^ 


experiment 


oo 


2.53 X lO-^mol/m^ 


[40] 


ki 


8.1 X lO^m^mol-is-i 


[41] 


k2 


6.7 X lO^m^mol-is-i 


[42] 


ks 


1.8 X lO-im^mol-is-i 


[28] 


rie 


2 


[33] 


A 


4.5 X 10"^ m^ 


experiment 



3.3. Governing Equations 



The governing equations for a chemical reaction network can be formulated by the law of mass 
action [4,21]. Coupling reactions in the enzyme layer with the one-dimensional-in-space diffusion. 
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described by Pick's second law, leads to the following equations of the reaction-diffusion type 

(0 < X < ai, t > 0): 



dt 

dercd 



-heoxSl + fc2ered"^ox,l (9) 
^iCoxSl - fc2ered"^ox,l (10) 



dt 

^ - D 
dt ^'^ dx 

"^ox,l j-^ "'ox,l 



^s.i^r^ - ^iCoxSi (11) 



= -PMox.I g^2 ' ~ ^2ered"^ox,l + /^3"^red,lOl (12) 



-PM.ed.l ' + fe2ercd"^ox,l - fcs^red.lOl (13) 

— = Doa,!^^ - fc3"^red,lOl (14) 

where x stands for space, t is time, eox(a;, t) and ei.ed(a;, t) correspond to concentrations of oxidized (Eqx) 
and reduced (Ered) enzyme, respectively; si(x, t) is the concentration of the substrate in the enzyme 
layer; mox,i(a;, t) and m,.cd,i{x, t) are the concentrations of the oxidized (Mqx) and reduced (M^d) forms 
of the mediator in the enzyme layer, oi is the concentration of oxygen in the enzyme layer and Ds^i, 
-Dmox,!' -^M.cd,!' -^02,1 are the diffusion coefficients of the corresponding substances defined by the 
subscript. In the definition of the concentrations and diffusion coefficients, here and later in this paper, 
the last numeric subscript label denotes the region of the model, particularly, 1 stands for the enzyme 
layer. The molecules of both enzyme forms, Eqx and Ercd, are considered as immobile, and therefore, 
there are no diffusion terms in the corresponding equations. 

Products Pi and P2 do not act as reactants in any reaction, so their concentrations are not used in any 
further calculations. Therefore, Equations (9)-(14) contain no equations for product Pi and P2. 

No enzyme molecules appear in other layers of the biosensor. Hence, only the mass transport by 
diffusion of species and the reaction (7) take place in other biosensor regions. The governing equations 
for these layers are represented as follows: (aj_i < a; < a^, t > 0, i = 2, 3, 4): 

^--^^ 

<9mox,i „ S^mox.j , , ,,,, 

= J-^Mo^,i—Q^ ^ ksmred^iOi (16) 

dmred,i d'^m^cd,i , 

—Qf- = DM,c^,i—g^2 hmred,iOi (IV) 



— = Do2,i-7^ - k3mred,iOi (18) 

where i = 2 corresponds to the PVA layer, 2 = 3 corresponds to the terylene membrane layer and i = A 
corresponds to the diffusion layer. 

3.4. Initial Conditions 



Let X = ao = 0 represents the electrode surface, while x = ai, x = a2, x = a-^ and x = a4 represent 
the boundaries between the adjacent regions, as described in Section 3.2 and shown in Figure 1. The 
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biosensor operation starts when the substrate and mediator appear in the bulk solution. This is used in 



the initial conditions {t = 0), 

e^-edix, 0) = 0, eox(x, 0) = Co, 0 < a; < Oi (19) 

Si{x, 0) = mox,j(a;, 0) = rriredA^, 0) = 0, a^.i < x < ai, i = 1,2,3 (20) 

S4(x,0) = mox,4(a;,0) =0, < x < (21) 

rriredAix, 0) = 0, as < X < 04 (22) 

54(04,0) = So, mox,4(a4,0) = mo (23) 

Oi{x,0)=oo, ai_i < X < fli, i = 1,2,3,4 (24) 

where Cq stands for the total concentration of the enzyme in the enzyme layer (eo = eoxix,t) + 



e^ed{x, t),Wx, t : X e {0,ai) ,t > 0), Sq stands for the substrate concentration, tjiq is the concentration 
of the oxidized form of the mediator and oq is the concentration of the oxygen in the bulk solution. 



3.5. Matching Conditions 



On the boundary between two adjacent regions having different diffusivities, matching conditions 
have to be defined (t > 0, z = 1, 2, 3) [26,27], 



D 



S,i 



dsj 
dx 



D: 



D 



x=a 

dm. 



ox, I 



dx 
drriredA 



S,i+l 



D 



dx 



Si{ai,t) = Si+i{ai,t) 



x=a 

(9m, 



Mox,j + l 



ox,i+l 



D 



Mox,i+l 



dx 
dm,.ed,i+i 



dx 



^ox,i{(J'i:t) — mox,i+l('3i, t) 



Do2,i 



doj 
dx 



doi+i 



dx 



Oi{ai,t) = Oi+i{ai,t) 



(25) 
(26) 
(27) 
(28) 



where i = 1 corresponds to the boundary between the enzyme layer and the PVA layer, i = 2 corresponds 
to the boundary between the PVA layer and the terylene membrane, whereas i = 3 corresponds to the 
boundary between the terylene membrane and the diffusion layer. 

These conditions mean that fluxes of the species through one region are equal to the corresponding 
fluxes entering the surface of the neighboring region. Concentrations of species in one region versus the 
neighboring region are assumed to be equal. 

3.6. Boundary Conditions 



In the bulk solution, the concentrations of species remain constant (t > 0), 

54(04, t) = So 
"^ox,4(a4,i) = mo 
m,edA{0'4,t) = 0 
04(04, t) = Oo 



(29) 
(30) 
(31) 
(32) 
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The reduced mediator, Mred, takes part in electrochemical Reaction (8) at the electrode surface 
(x = 0). The rate of this reaction is considered so high that the concentration of Mred at the electrode 
surface is permanently reduced to zero (t > 0), 



"^red,l(0)^) = 0 



(33) 



Since electrochemical Reaction (8) produces as much Mqx as it consumes Mred, the flux of Mqx on 
the electrode surface is equal to the flux of Mrcd, but in the opposite direction. This relation is expressed 
by the following boundary condition (t > 0): 



dm 



red,l 



dx 



--DMox,1 



<9mox,i 



x=0 



dx 



(34) 



x=0 



The substrate, S, and oxygen, O2, are electrode-inactive species; thus, their fluxes on the electrode 
surface are equal to zero (t > 0), 



dsi 
dx 
doi 
dx 



x=0 



x=0 



0 
0 



(35) 
(36) 



3. 7. Biosensor Response 



The measured current is usually assumed as the response of an amperometric biosensor in physical 
experiments. The biosensor current, i{t), at time t was expressed explicitly from Faraday's and 
Pick's laws, 

dm^^Ai 



i{t) = ArieFD 



M„H,1 



dx 



(37) 



x=0 



where i{t) is the faradaic current generated by electrochemical Reaction (8), A is the electrode surface. 
Tie is the number of electrons involved in a charge transfer at the electrode surface and F is the Faraday 
constant, F = 96, 486 C/mol. 

We assume that the system approaches a steady state as t — )• 00, 



(-St 



lim i(t) 



(38) 



where i^t is the the steady-state biosensor current. 

Let us introduce Bq^, which shows the biosensor response sensitivity to oxygen. 



02 



istjop = 0) - istjoo = Co J 
ist{oo = 0) 



(39) 



where ist{oo = 0) is the biosensor response in anaerobic conditions, ist{oo = cq^) is the biosensor 
response in aerobic conditions and cq^ = 2.53 x 10^'' mol/m^ is the oxygen concentration in water [40]. 
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4. Numerical Simulation of Biosensor Action 



4.1. Simulating the Biosensor Operation 



The exact analytical solution for the problem (9)-(36) is not known. Therefore, the problem was 
solved numerically, using the finite difference technique [43,44]. An implicit finite difference scheme 
was built on a uniform discrete grid with 50 points in the space direction for each modeled layer 
corresponding to a certain time moment. The simulator has been programmed by the authors in the 
C++ programming language [45]. 

In the numerical simulation, the biosensor response time was assumed as the time when the change 
of the biosensor current remains very small during a relatively long term. A special dimensionless decay 
rate, e, was used. 



mm 

j(i)>0 



t : 



t 

W) 



di(t) 



< £ 



where is the biosensor response time. The decay rate value e 



i{tr) ~ ist (40) 
= 10^^ was used in the calculations. 



Figure 2. Experimental (1, 2, 3) and computational (4, 5, 6) calibration curves at 
different concentrations of the mediator (mo): 0.2mol/m'^ (1, 4), 0.05mol/m^ (2, 5), 
0.005 mol/m^ (3, 6). 




Sq (mol m"^ ) 



4.2. Model Validation 



The numerical solution of the model (9)-(36) was compared with the experimental data. The results 
are depicted in Figures 2 and 3. 

As one can observe from Figure 2, the simulated calibration curve at mo = 0.05mol/m^ fits the 
experimental data well. At higher and lower mediator concentrations, the modeled data do not fit the 
experimental results so well. This may be explained by some processes that are not accounted for in the 
mathematical model. These processes may include reverse reactions, the instability of the compounds or 
enzyme degradation. Furthermore, some reaction rate constants, ki, k2 or k^, may exhibit dependence 
on the mediator concentration. 

The dynamics of biosensor responses during the experiments was compared with the results of the 
computational experiments. The data are depicted in Figure 3. 
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Figure 3. The dynamics of the biosensor current at different values of the substrate 
concentration, sq: 0.49mol/m3 (1), 0.99mol/m3 (2), l.QQmol/m^ (3), 4.98mol/m3 (4), 
9.9mol/m'^ (5), 19.6mol/m'^ (6); (a) experiment; (b) simulation. 
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As one can observe from Figure 3, the dynamics of the simulated biosensor response is in good 
agreement with the experimental data. At lower concentrations of the substrate (sq = 0.49mol/m'^ and 
So = 0.99mol/m^), the current difference at the steady state is 6.1% and 7.5%, respectively. At the 
intermediate concentrations (sq = 1.99mol/m'^ and sq = 4.98mol/m^), the current difference at the 
steady state is slightly higher (9.7% and 18%, respectively). At higher concentrations of the substrate 
(so = 9.9 moXlvD? and sq = 19.6 moVm^), the differences in the steady-state current are 4.6% and 0.58%, 
respectively. These results show that the simulated data more accurately coincide with the experimental 
data at relatively low and high substrate concentrations. 

4.3. Concentration Profiles at the Steady State 

Concentration profiles may be a useful source of information about the processes that take place 
during the biosensor operation. The biosensor action was simulated at the parameter values as specified 
in Table 1. The concentration profiles at the steady state are depicted in Figure 4. 

One can observe in Figure 4 that the concentration of oxygen is almost the same through all the 
investigated region. The oxygen concentration at the electrode surface is only 3.6% lower than in the 
bulk solution. This indicates that the rate of Reaction (7) is comparatively slow in comparison with the 
diffusion of this gas. 

The curve representing the concentration of Mred shows that the concentration at the electrode surface 
is reduced to zero, due to the fast electrochemical Reaction (8). The same electrochemical reaction 
causes an increase in the concentration of Mqx at the electrode surface, where it has practically the 
same concentration as in the bulk. In the middle region, the concentration of M^ed is higher than the 
concentration of Mqx- 

The concentration of the substrate at the electrode surface is equal to 3.93mol/m^, which is 21% 
less than in the bulk solution. This significant reduction in the concentration shows that the rate 
of Reaction (5) is comparatively high in comparison with the diffusion of the substrate from the 
bulk solution. 
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Figure 4. The concentration profiles of S (1), E^ed (2), Eox (3), M,ed (4), (5) and O2 (6) 
at the steady state. Vertical axes correspond to the concentration of species, horizontal axes 
correspond to the distance from an electrode surface. The vertical lines mark the boundaries 
of the biosensor layers. Concentrations in graph (a) are provided in logarithmic scale. In 
graph (b) concentrations of four species are depicted on a linear scale. 
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One can observe that the concentrations of Ej-cd and Eqx correlate with the concentrations of M^ed 
and Mox, respectively. While in the middle regions, the E^d concentration is higher than Eqx, it is the 
opposite at the electrode surface. The concentration shift at the electrode surface may be explained by 
the fact that Ercd is the reactant in Reaction (6), the rate of which is significantly increased by the increase 
in the Mqx concentration, which is the second reactant of Reaction (6). 



4.4. The Biosensor Response Sensitivity to Oxygen 

In order to assess the biosensor response sensitivity to oxygen, the biosensor operation was simulated 
in both aerobic and anaerobic conditions. Bq^ was used as a measure of the sensitivity. 

In order to investigate what properties a compound should have to be a suitable mediator, the 
biosensor response sensitivity, -B02 » dependence on ^3 was investigated at different values of the substrate 
concentration, sq. The other biosensor parameters were equal to the values indicated in Table 1. The 
results are depicted in Figure 5. 

The values of Bq^ depicted in Figure 5 indicate that the design of the biosensor is very successful. The 
simulated biosensor response sensitivity to oxygen shows that the oxygen influence is about 6-7 orders 
of magnitude smaller than the biosensor response itself. From the viewpoint of biosensor sensitivity 
to oxygen, the PMS is a very successful mediator, as the value of Bq^ is very low at the value 
corresponding to PMS. These results are in accord with physical experiments, which did not observe a 
difference in readings while conducting the same experiment in both anaerobic and aerobic conditions. 

Even though the biosensor response sensitivity to oxygen is small, it is dependent on the value of 
^3 and substrate concentration sq. Two main tendencies may be observed from the curves depicted in 
Figure 5: the oxygen has a bigger influence when is bigger and when the substrate concentration, sq, is 
smaller. Bigger values of reaction rate constant k^ results in higher rates of Reaction (7). Therefore, this 
directly influences the value of Bq^ . However, if the concentration of the substrate is high, the change 
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in ^3 only slightly alters the values of Bq^. At lower substrate concentrations, the shift towards higher 
biosensor response sensitivity to oxygen is more significant if k^, increases. 

Figure 5. The biosensor response sensitivity, Bq^, to oxygen dependence on reaction rate 
constant of Reaction (7), k^, at different values of substrate concentration sq: 0.498 mol/m^ 
(1), 4.98 moMvD? (2) and 49.8 mol/m^ (3). The vertical line marks the reaction rate constant, 
/c3, corresponding to N-methylphenazonium methyl sulfate (PMS). 




10"^ 10"^ 10"' lo" lo' 10^ 10^ lO"^ 
^3 (m'' mol"' s"') 



The rate of Reaction (6) is also dependent on the properties of a mediator; particularly, it is directly 
dependent on the reaction rate constant, k2. The biosensor response sensitivity, Bq,^, to oxygen 
dependence on k2 was investigated at different values of the substrate concentration, sq. The results 
are depicted in Figure 6. 

Figure 6. The biosensor response sensitivity, Bq^, to oxygen dependence on reaction rate 
constant of Reaction (6), k2, at different values of substrate concentration sq: 0.498 mol/m^ 
(1), 4.98 moMw? (2) and 49.8 mol/m^ (3). The vertical line marks the reaction rate constant, 
/c2, corresponding to PMS. 
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As is evident from Figure 6, the biosensor response sensitivity to oxygen is even less dependent on 
the reaction rate constant, k2. The main tendency that may be observed from Figure 6 is that biosensor 
sensitivity Bq.^ is lower at higher values of reaction rate constant ^2. The sensitivity of the biosensor to 
oxygen at /c2 value corresponding to PMS is one of the lowest in the range of investigated sensitivities. 
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However, besides the main tendency, one may observe that at lower substrate concentrations 
(so = 0.498 mol/m'^), the dependency is non-monotonous. The biosensor response sensitivity to oxygen 
is the lowest when the value of reaction rate constant ^2 is equal to approximately 5 x 10"^ m'^ mol~^ s^^. 
At lower and higher values of /c2, the biosensor shows a higher sensitivity to oxygen. 

The impact of the diffusion layer thickness, ^4, on the biosensor response sensitivity to oxygen was 
also investigated. Simulation results are depicted in Figure 7. 

Figure 7. The biosensor response sensitivity, Bq^, to oxygen dependence on the 
thickness, ^4, of the diffusion layer at different values of the substrate concentration, 
So:0.498mol/m3 (1), 4.98mol/m3 (2) and 49.8mol/m3 (3). The vertical line marks the 
estimated diffusion layer thickness during the physical experiment. 
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The curves in Figure 7 show that the impact of diffusion layer thickness on the biosensor sensitivity 
to oxygen is small. However, the tendency that a thicker diffusion layer results in a higher oxygen 
influence may be observed. This may be explained by the fact that the diffusion coefficient of oxygen is 
the highest among the compounds that are reactants in reaction scheme (5)-(8) (see Table 1). The thicker 
the diffusion layer, the faster that oxygen reaches the biosensor compared with the other compounds. 

5. Conclusions 

The mathematical model (9)-(36) describes the processes that take place during the biosensor 
operation sufficiently well. The model may be used as a tool for assessing the need to conduct the 
measuring experiment with a biosensor in anaerobic conditions. The model may be also used to 
investigate other properties and behaviors of a biosensor. 

The comparison of experimental and simulated results showed that the model most accurately 
describes the biosensor operation at moderate mediator concentrations (mo = 0.05 mol/m'^). At lower 
and higher concentrations of the mediator, the accuracy is lower. 

The investigation of simulated concentration profiles at the steady state showed that the consumption 
of oxygen is very low. This shows that the rate of reaction consuming molecular oxygen is low 
and should not significantly influence the biosensor response when an experiment is carried out in 
aerobic conditions. 

The investigation of the biosensor response at different parameter values showed that overall, the 
biosensor response sensitivity to oxygen is very small. The change in response in anaerobic conditions 
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is six to seven orders of magnitude smaller than the biosensor response itself. It was shown that the 
mediator was very successfully chosen from the viewpoint of biosensor response sensitivity to oxygen. 

It was also shown that a thicker diffusion layer slightly increases the biosensor response sensitivity to 
oxygen. This is consistent with the fact that the diffusion coefficient of oxygen is the highest among the 
compounds participating in chemical reactions during biosensor operation. 

Simulated results and conclusions inferred from them should be used with caution, as the biosensor 
model poorly fits the experimental data with some mediator concentrations. This may indicate that some 
processes are not accounted for in the biosensor model or that some reaction rate constants used in the 
model exhibit dependence on the mediator concentration. 
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